A^— from Renormalization Group Optimized Perturbation 
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A recent extension of a variationally optimized perturbation, combined with renormalization group 
properties in a straightforward way, can provide approximations to nonperturbative quantities such 
as the chiral symmetry breaking order parameters typically. We apply this to evaluate, up to third 
order in this modified perturbation, the ratio F n / A-jjg, where F n is the pion decay constant and 

A-jjg the basic QCD scale in the MS scheme. Using experimental F„ input value we obtain A^L. 2 ~ 

255+1° MeV > where quoted errors are estimates of theoretical uncertainties of the method. This 
compares reasonably well with some recent lattice simulation results. We briefly discuss prospects 
(and obstacles) for extrapolation to as{.n) at perturbative \x values. 



In the chiral symmetric, massless quarks limit, the 
strong coupling as(fJ-) at some reference scale \i is the 
only QCD parameter. Equivalently the Rcnormalization- 
Group (RG) invariant scale 
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in a specified renormalization scheme, is the fundamental 
QCD scale. In (1) /3o, /3% are (scheme-independent) one- 
and two-loop RG beta function coefficients, and ellipsis 
denote higher orders scheme-dependent RG corrections 
as will be specified below. As indicated A— — also de- 
pends on the number of active quark flavors n/, with non- 
trivial (perturbative) matching relations at the quark 
mass thresholds (see e.g. the QCD chapter in ([1]) for 
a recent review and original references) . as has been ex- 
tracted from many different observables confronted with 
theoretical predictions, and its present World average is 
impressively accurate [1]: as{mz) — .118±.001 (though 
there are long-standing tensions with values from struc- 
ture functions [2]: as{mz) — .114 ± .002). In any case, 
it is still of great interest to estimate A-p-^ from other 
observables, and other theoretical methods, specially to 
access the deep infrared, nonperturbative QCD regime 
for rif — 2(3), where a perturbative extrapolation from 
otsip^z) is unreliable. Indeed for several years determi- 
nations of A— — for nt < 2(3) from Lattice simulations 

M S — 

have been the subject of much activities. 

In this letter we explore a different route to estimate 
such quantities, more rooted in perturbation theory, and 
where the dynamically broken chiral symmetry due to the 
light u, d (and s) quarks plays a crucial role. The main 
order parameters of chiral symmetry breaking, namely 
the chiral quark condensate (qq) and pion decay con- 
stant Ftj-, should be entirely determined by the unique 
scale A-pg in the strict chiral limit. Well-established ar- 
guments usually consider hopeless to calculate the above 
order parameters from QCD first principle, except on 
the lattice. First, most obviously because of the men- 
tioned nonperturbative regime at the relevant scale close 
to A-p-g, implying a priori large as values invalidating 



reliable perturbative expansions. Second, standard per- 
turbative series of those quantities at arbitrary orders 
are anyhow proportional to the quark masses m q (up to 
powers of lnm g ), so trivially vanish in the strict chiral 
limit m q — > 0. Moreover, general arguments, related to 
the problem of resumming presumed factorially divergent 
perturbative series at large orders [3], seem to further in- 
validate any perturbative approach to calculate the order 
parameters. We will see how to circumvent at least the 
first two arguments above, by a peculiar modification of 
the ordinary perturbative expansion in as, with possi- 
ble improvements of the convergence and resummation 
properties of such a modified series. In this letter we con- 
centrate on determining F^/ Ajjg at successive orders of 
this modified perturbation, thus extracting Aj^g values 
from the pion decay constant value F v . 
Our method has been recently applied [4] to the D = 2 
Gross Neveu (GN) 0(2N) model[5], which shares many 
properties with D = 4 QCD: it is asymptotically free, 
has a (discrete) chiral symmetry for m = 0, dynamically 
broken with a fermion mass gap. The exact mass gap 
is known for arbitrary N, from Thermodynamic Bcthc 
Ansatz[6], allowing accurate tests of our method. Using 
only the two-loop ordinary perturbative information, we 
obtained approximations to the exact mass gap at the 
percent or less level [4], for any N values. 

The basic framework [7] is to introduce an unphysical 
parameter < 5 < 1, interpolating between Cf ree and 
Winter action , such that the relevant fermion (quark) mass 
m q becomes an arbitrary "variational" parameter: 

C QC D(m q ,as) -> C Q c D (m(l - S) a ,a s 8) (2) 

where £qcd((xs) stands for the standard complete QCD 
Lagrangian, and m q originally is a current quark mass 
relevant for chiral symmetry breaking. In the following 
we shall mainly consider two quark flavors u, d and the 
corresponding SU(2) L x SU(2) R — > SU{2) V chiral sym- 
metry breakdown, with m = (m u + ma)/ 2 neglecting as 
usual the m u — rrid difference. The extra parameter a in 
(2) reflects the large freedom in the modified interpolat- 
ing Lagrangian, and will allow imposing further physical 
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(or technical) constraints, as we shall specify later. 
The procedure is fully consistent with renormalizability 
and gauge invariance, provided that the above redefini- 
tion of the QCD coupling as —> Sag is performed consis- 
tently for all interaction terms appropriate for gauge in- 
variance and renormalizability. Working with the above 
Lagrangian is perturbatively equivalent to taking any 
standard renormalized series in g = 4iras for a physical 
quantity, re-expanded in powers of 5 after substitution: 



m ->■ m (1 - 6) a , g^5g. 



(3) 



One takes afterwards the 5 — > 1 limit to recover the 
original massless theory. This expansion gives, how- 
ever, a remnant m- dependence at any finite <5 fe -order, 
and m can be fixed conveniently by an optimization 
(OPT) prescription [8]. The convergence of such a pro- 
cedure, which may be viewed particular case of 
"order-dependent mapping" [9], has been proven[10] for 
the D = 1 A0 4 oscillator model. In renormalizable D > 1 
models, the situation is more involved and it is difficult to 
make statements on the possible convergence properties 
(see however [11] for a particular case). But at least the 
method allows to obtain approximations to nonperturba- 
tive quantities beyond the mean field approximation in 
various models, which (empirically) appear to converge 
rather quickly at the first few perturbative orders. 

Previous attempts to use this approach in QCD gave 
rough estimates of the order parameters (dynamical 
"mass gap", F„, (qq))[12}. But it involved a cum- 
bersome manner of incorporating renormalization group 
(RG) properties within such modified perturbative series, 
difficult to generalize beyond the first or second RG order 
and to other physical quantities defined by their pertur- 
bative series. Our new proposal introduces in contrast 
a much simpler marriage of OPT and RG properties [4]: 
consider an ordinary perturbative expansion for a phys- 
ical quantity P(m,g), after applying (3) and expanding 
in 5 at order k. In addition to the OPT equation: 



/-P( fc >(m, ff ,5 = l)| r 
am 



0, 



(4) 



we require the (<5-modified) series to satisfies a standard 
RG equation: 
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where the usual RG operator 

d d n . . d , . d 



(5) 



(6) 



gives zero to 0(g k+1 ) when applied to RG-invariant 
quantities. (NB our normalization is (3{g) = dg/dln^ = 
-2b g 2 -2b 1 g 3 + --- and 7m (.g) = j a g + jig 2 + - • •■ The b t 
and -fi known up to 4-loop are given in [13]). Note that, 



combined with Eq. (4) , the RG equation takes a reduced 
form: 



dg 



P {k) {m,g,5 =1) = 



(7) 



and Eqs. (7), (4) completely fix [for given a values in 
Eq. (3)] optimized values m = m and g = g. 

We shall now illustrate the method concretely on a 
well-defined perturbative series relevant for the pion de- 
cay constant F v . One very convenient definition of F^ is 
via the axial current correlator, known at present up to 
4-loop orders [14, 15]. More precisely: 

i(0|T4(p)^(0)|0> = S^g^F 2 + 0(p, Pv ) (8) 

where the axial current is A % = 97^,75 ^ q, and in this 
normalization F n ~ 92.3 MeV [1]. 

Our starting point is thus the perturbative expansion of 
(8) in the MS scheme: 

Ffoert) = 3f£ [div(e, a s )-L+ ff (8i 2 + \L + |) 



+(ff ) 2 [/ 30 i 3 + fsiL 2 + f 32 L + f 33 ] + 0{a%) 



(9) 



where L = In : 



f _ 304 32^, f _ 136 , 32^, 

M , J30 - — - —n f , j 31 - — 3- + T n f , 
and fs2 and the non-RG coefficient J33 have more lengthy 
expressions easily extracted from related calculations in 
[14] valid for arbitrary numbers of quark flavors. Re- 
cently even the 0(o? s ) coefficients fu,i — 0, • • • ,4 were 
obtained [15], that we also use in our analysis [16]. 

There is however one subtlety at this stage: as is 
well-known, at this level the calculation e.g. in dimen- 
sional regularization of (8) actually still contains diver- 
gent terms after mass and coupling renormalization in 
MS scheme, formally indicated as div(e, as) in Eq. (9). 
This simply reflects the extra additive renormalization 
needed for such a composite operator. But to obtain a 
RG-invariant finite expression from (9) the subtraction of 
those divergences should be performed consistently with 
RG properties. Now to fix this subtraction at order k 
needs knowledge of the coefficient of the L term (equiva- 
lently the coefficient of 1 /e in dimensional regularization) 
at order k + 1. We define the needed subtraction as a per- 
turbative series: 



sub(g,m) = m 2 ^s l g l 



-1 



(10) 



i>0 



with coefficients determined order by order by 



li— [sub(<7, m)} = Remnant(<7, m) (11) 

where the remnant part is obtained by applying the RG 
operator Eq. (6) to the finite part of (9), as the latter is 
not separately RG-invariant. Thus the (finite) quantity 
F 2 (pert) (finite) — sub(g,m) is RG-invariant at a given 
order. Note that (10) does not contain any L terms 
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and necessarily starts with a So/g term to be consis- 
tent with RG invariance properties. (This reflects the 
fact that the one-loop contribution in Eq. (9) is of order 
g°). Completely equivalent results are obtained [12] more 
formally by working with bare expressions and establish- 
ing the required RG properties. We obtain for instance 

S = 4^ (bo -70) ' Sl = leffi-Q+an,) ' &nd hi § her 0rder S * 

have more lengthy expressions not given here. 

We thus apply to (9)-(10) the procedure (3) and ex- 
pand at order 5 k , then solving OPT and RG Eqs.(4), 
(7). Before coming to numerical results, some important 
remarks are in order. First, Eqs.(4), (7) being polyno- 
mial in [L,g), one serious drawback is that at increasing 
5-orders there are (too) many solutions, most being com- 
plex (complex conjugate in fact since all coefficients of 
(4), (7) are real). Now an important selection comes 
about if one imposes an additional constraint that the 
solutions should obey asymptotically the standard per- 
turbativc RG behavior for g — »■ 0: 



TABLE I. Combined OPT+RG results at successive 5-order 



g(fj, > to) 



(26 ln^) 

TO 



(12) 



This is a very natural requirement, otherwise optimal 
solutions do not match standard perturbative behavior. 
An important related remark is that, after OPT, the 
optimized mass to is consistently O(Ajjg) (rather than 
to ~ 0). Thus to plays the role of a mass gap, such 
that the OPT-modihcd expansion F v oc to is a pertur- 
bation around a 'Born-level' value of O(A-p-g) in con- 
trast with the original standard perturbative expansion, 
F n (m — > 0) — > 0. In fact, at 5 fe -order, Eq. (7) is a polyno- 
mial of order k + 1 in L, thus exactly solvable up to third 
order, with full analytical control of the different solu- 
tions. However, unlike the GN mass gap case [4], to have 
at least one of the RG OPT F w solutions behaving as 
(12) at any <5 fc -orders, requires a critical value of the pa- 
rameter a in the interpolation (3), namely a = 7o/(26o). 
This connection with RG anomalous dimensions is not 
too surprising: similarly in other theories, e.g. <3> 4 in 
D = 3, specific a values occur, consistent with RG crit- 
ical exponent properties, as emphasized in [17], at the 
same time matching real optimal solutions [18]. 
We thus fix a = 7o/(2&o) to determine solutions at suc- 
cessive (5-orders. Most solutions exhibit very odd depen- 
dence in g incompatible with (12), as also observed very 
similarly in the GN model case [4] . This perturbative RG 
criteria appears to give unique solutions (at least up to 
third order here considered), given in Table I. 

But those unique well-behaving RG solutions remain 
complex (conjugates) for F n . It can always be that other 
models, other physical quantities [4, 18] (or very different 
rif values) give real solutions, but not for F n at k < 3 
orders here explored. Since this is unphysical, we can 
only expect acceptable solutions of behavior (12) to have 
at least Re(g) > and lm(F n ) <C Re(F 7r ), the imaginary 
part indicating an intrinsic theoretical uncertainty of the 



(5-order k 


K (m,g) 
A PA 
MS 


L 


a s 


1 


0.372 ±0.16i 


-0.45 ±0. Hi 


1.01 ±0.08i 


2 


0.353 ± 0.03i 


-0.52=F0.69i 


0.73 ± 0.02i 


3 (s 4 , /44 = 0) 


0.351 ± 0.08i 


-0.13=F0.04i 


0.61 ±0.33i 


3 (s 4 = PA[1,2]) 


0.341 ± 0.07i 


-0.23=F0.04i 


0.59±0.31i 



method, as will be specified below. 

To compare our results with other (principally Lattice) 
calculations, one should be careful to use the same con- 
ventions for Aj^g. We mainly use a convenient (Pade 
Approximant) 3- loop form, cf. [19]: 



MS ' 



' 2b s 



, 1 ^ ^ bo bj 



(13) 



We also compare with a more standard 4-loop perturba- 
tive form [1], with 63 ^ 0, which gives a systematic ~ 2 
% lower A-p-g values for our optimal as values. 

Comparing second and first 5-orders in Table I, one 
observes that the solution has a much smaller imaginary 
part, and also Re as decreases to reasonably perturbative 
values as the (5-order increases. At third order, the g 3 S4 
term in (10) needs knowledge of the presently unknown 
5-loop coefficient of L. We have thus estimated S4 either 
with a Pade Approximant PA[1, 2] from lower orders, or 
alternatively simply ignoring S4,/44 = 0, retaining only 
4-loop RG ln p (TO.//n) coefficients. The difference between 
those two choices in Table I gives one estimate of higher 
order uncertainties. We also incorporate additional theo- 
retical uncertainties by solving Eq. (7) truncated to lower 
g orders, or neglecting 63, etc (since RG- invariance is 
only required up to 0(g k+1 ) terms at order k). Op- 
timal RG solutions are remarkably stable with respect 
to such approximations on 4-loop order and RG trun- 
cations, with at most ~ 2-3% differences on . In 
addition, as above mentioned we take into account a more 
intrinsic error: given the (unphysical) imaginary parts of 
the solutions, we empirically take the range spanned by 
Re(F^(g,L)) — F n (Ke(g), Re(Z)), as this tends to maxi- 
mize the uncertainty for increasing Im(g,L). This gives 
only about a 1-2% variation on at 0(S 2 ); but a 
larger ~ 10% one at 0(5 3 ) due to the larger imaginary 
part of the solutions, perhaps an artefact of the unknown 
exact S4 coefficient at this order. Since optimal solutions 
are complex conjugates, another estimate could be sim- 
ply to compare their real parts with their modulus, which 
gives a much more moderate difference (2% at 0(S 3 )). 
Clearly the occurrence of complex solutions is our main 
source of theoretical uncertainties, so there is potentially 
room for improvements, e.g. from other more general 
prescriptions [18]. However, we prefer to keep a conser- 
vative estimate of theoretical errors at this stage. 
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Finally we can subtract out the explicit chiral symme- 
try breaking effects from small m u , nid ^= 0: it is in princi- 
ple possible to incorporate those effects within the varia- 
tional framework [12], but for the time being we shall sim- 
ply rely on other known results. Defining F as usual as 
the F n value in the strict chiral limit m u , raj —> 0, Lattice 
simulations recently obtained [20]: -^f- ~ 1.073 ± 0.015, 
that we accordingly take into account in the final A^T 2 
numerical value. With all theoretical uncertainties (lin- 
early) combined we obtain: 



A^£r 2 ~ 255 ± 15+ 25 MeV . 

MS 



(14) 



The central value corresponds to Re F 2 (g, L), the first er- 
rors encompass both higher order and F^/F above men- 
tioned uncertainties, while the upper bound corresponds 
to^ 2 (Rc(5),Re(L)). 

One may compare this with three main classes of lattice 
calculations based on very different methods. First in 
the Schrodinger functional scheme[21], A^T 2 = 245 ± 
16(stat) ± 16(syst) MeV. Next for Wilson fermions [19]: 
A^T 2 = 261 ±17 ±26 MeV. Finally for twisted fermions 
(including nonperturbative power corrections in analysis) 
[22]: A^T 2 = 330 ± 23 ± 22_ 3 3 MeV. Those differences 
are presumably related to different dynamical quark mass 
values in different Lattice calculations, and also differ- 
ent chiral extrapolation methods (see e.g. the discussion 
in [22]). 

Finally we could in principle extrapolate to as(p) at 
high (perturbative) scale /x. A main obstacle however is 
going from nf — 2 to nf — 3, crossing the strange quark 
m s threshold in the deep infrared regime, where one can- 
not use standard perturbative extrapolation. But (9) be- 
ing known for arbitrary n/, we can calculate similarly 
Fn/Ajt^ 3 . The details arc skipped for elsewhere, but the 

outcome is a mild variation, with F V /A^LZ > F^/A^T 
by only a few percent. However, the final value of 
A^T is much dependent also on the ratio F^/Fo, where 
Fo = F^{m u , nid, m s — > 0). The amount of explicit chiral 
symmetry breaking from m s ^ is clearly more impor- 
tant than in the nj — 2 case, and indeed still subject to 
intense debates, with still large uncertainties even from 
Lattice results [20]. Moreover, since our RG-improved 
OPT modifies perturbative expansions, it should also be 
used consistently to extrapolate to higher scales, which 
can differ substantially from a standard perturbative ex- 
trapolation. For both those reasons we refrain from giv- 
ing a precise prediction of asijnz) at this stage, mainly 
due to the large uncertainties involved in subtracting out 
explicit chiral symmetry breaking from m s . We plan in 
the next future to implement those effects directly within 
the OPT framework. 



In conclusion, a straightforward implementation of RG 
properties within a variationally optimized perturbation 
was proposed, using only perturbative information at the 
first few orders. In QCD, calculations at first (5-order give 
already very reasonable approximations to F n / Ajj^, and 
second and (approximate) third order results exhibit a 
remarkable stability. These results compare reasonably 
well with recent lattice calculations of Aj^, though the 
best with those in [19]. We conjecture that the remain- 
ing discrepancies (with other lattice results, and possibly 
with the World average «s(mz) values) could be due to 
the interplay with explicit quark mass effects, which are 
in principle implementable in our framework. The precise 
extrapolation to as{mz) is however beyond the present 
scope and postponed for a future work. 
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